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Abstract 

The paper gives an introduction to the gravitational radiation theory of isolated sources and 
to the propagation properties of light rays in radiative gravitational fields. It presents a theo- 
retical study of the generation, propagation, back-reaction, and detection of gravitational waves 
from astrophysical sources. After reviewing the various quadrupole-moment laws for gravita- 
tional radiation in the Newtonian approximation, we show how to incorporate post -Newtonian 
corrections into the source multipolc moments, the radiative multipole moments at infinity, and 
the back-reaction potentials. We further treat the light propagation in the linearized gravita- 
tional field outside a gravitational wave emitting source. The effects of time delay, bending of 
light, and moving source frequency shift are presented in terms of the gravitational lens potential. 
Time delay results are applied in the description of the procedure of the detection of gravitational 
waves. 

Pacs Numbers : 04.25.-g, 04.25.Nx 



1 Introduction 

It was only in the late fifties of the twentieth century that by the work of Hermann Bondi and 
Joseph Weber gravitational radiation entered the domain of physics. Before that time gravitational 
radiation was not considered to be of observational relevance and the gravitational radiation theory 
was not developed very deeply. 

The supposed detection of gravitational radiation by J. Weber in the late sixties triggered strong 
and still on-going efforts both in the building of gravitational wave detectors and in the elaboration of 
the gravitational radiation theory, including investigations of the most reliable sources of detectable 
gravitational waves, calculations of wave forms, and analysis of data from detectors (cf. IQ). It 
turned out that coalescing neutron stars and/or stellar-mass black holes, together with gravitation- 
ally collapsing objects (type II supernovae), are the most relevant sources for detectable gravitational 
waves on Earth because they are strong and fit well to the frequency band of the Earth-based de- 
tectors which ranges from 10 Hz to 10 kHz. The strength of these sources is such high that several 
detection events per year might be expected in future fully developed detectors. The most sensitive 
Earth-based detectors will go into operation in the first few years of the new millenium. These 
are the laser- interferometric detectors in Germany, GEO600, built by a German/British consor- 
tium, in Italy, VIRGO, built by a Italian/French consortium, the two LIGO detectors in the United 
States, and the TAMA300 detector in Japan. There are several bar detectors already operating on 
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Earth (ALLEGRO in the United States, AURIGA and NAUTILUS in Italy, EXPLORER at CERN, 
NIOBE in Austraha). These detectors are being permanently upgraded and will be supplementing 
the measurements of the interferometric detectors later. For the measurement of gravitational waves 
in the frequency range between 0.1 Hz and 0.1 mHz the space-borne laser-interferometric detector 
LISA is devised which is expected to be flown around 2010 by NASA/ESA. The astrophysical sources 
of the gravitational waves to be detected by this detector are a variety of orbiting stars (interacting 
white dwarf binaries, compact binaries), orbiting massive black holes, as well as the formation and 
coalescence of supermassive black holes. Stochastic gravitational waves from the early universe are 
expected to exist in the whole measurable frequency range from lO'' Hz down to 10"^*^ Hz. The tools 
to possibly measure the primordial waves are the mentioned Earth-based and space-borne detectors, 
Doppler tracking, pulsar timing, very long baseline interferometry, as well as the cosmic microwave 
background. 

On the theoretical side there are essentially two approaches which permit to investigate the prop- 
erties of and to make predictions about gravitational waves from various sources. The first approach, 
that we can qualify as "exact" , stays within the exact theory, solving or establishing theorems about 
the complete non-linear Einstein field equations. Within this approach one can distinguish the work 
dealing with exact solutions of the field equations in the form of plane gravitational waves, and 
especially colliding plane waves. Since the waves are never planar in nature, this work is not very 
relevant to real astrophysics, but its academic interest is important in that it permits notably the 
study of the appearance of singularities triggered by collisions of waves. Also within the exact 
approach, but more important for applications in astrophysics, is all the work concerned with the 
study of the asymptotic structure of the gravitational field of isolated radiating systems. The work 
on asymptotics started with the papers of Bondi et al. and Penrose [^. The second approach is 
much more general, in the sense that it is not restricted to any particular symmetry of the system, 
nor it is applicable only in the far region of the system. However, the drawback of this approach 
is that it is only approximate and essentially looks for the solutions of the Einstein field equations 
in the form of formal expansions when c — > oo (post-Newtonian approximation). This approximate 
post-Newtonian method can be applied to the study of all theoretical aspects of gravitational radi- 
ation: the equations of motion of the source including the gravitational radiation reaction (works 
of Einstein, Infeld, and Hoffmann Chandrasekhar and Esposito ||5|, Burke and Thorne |^, |^, ^, 
Ehlers [^, Papapetrou and Linet |jlO|, Damour and Deruelle |11, Il2|, Schafer [13|, Kopejkin [1^]); 
the structure of the radiation field (work of Bonnor |15], Thorne [|l6[, Blanchet and Damour [17]), 
and, more recently, accurate post-Newtonian wave generation formalisms |18, 19, |2^, p^ ]. 

To the lowest, Newtonian order, the wave-generation formalism is called the quadrupole formal- 
ism, because as a consequence of the equality of the inertial and gravitational mass of all bodies the 
dominant radiating moment of any system is the (mass-type) quadrupole, which simply is at this 
approximation the standard Newtonian quadrupole moment. We are very much confident in using 
the theoretical framework of the post-Newtonian approximation because, marvellously enough, the 
framework of the Newtonian, quadrupole formalism has been checked by astronomical observations. 

In fact, there are two observational tests of the validity of the quadrupole formalism. The first 
test concerns the famous Hulse-Taylor binary pulsar whose decrease of the orbital period by 
gravitational radiation is predicted from the quadrupole formula to be ||2^, 24, |2^ 

■ _ 192^ /27rGy/=^ MpM, 1 + je' + je' . . 
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where Mp and are the pulsar and companion masses, e is the orbit eccentricity, and G, c are the 
universal gravitational constant and the speed of light. Numerically, one finds Pb = —2.4 x 10~^^ 



sec/sec, in excellent agreement (0.35% precision) with the observations by Taylor et al. 25]. The 
second test concerns the so-called cataclysmic variables. There we have binary systems in which 
a star filling its Roche lobe (the "secondary" with mass M2) transfers mass onto a more massive 
white dwarf (the "primary" with mass Mi > M2). From the formula for the angular momentum in 
Newtonian theory J = GMiM2(a/GM)^/^ (where M = Mi + M2), we deduce the secular evolution 



of the orbital semi-major radius a (whatever may be the mechanism for the variation of J), 



^ ' - ^ \ (2) 
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where M2 is the rate at which the secondary transfers mass to the primary (M2 < 0). Since 
Ml > M2, the mass transfer tends to increase the radius a of the orbit, hence to increase the radius 
of the secondary's Roche lobe, and, thus, to stop the mass transfer. Therefore a long lived mass 
transfer is possible only if the system looses angular momentum to compensate for the increase of a. 
For cataclismic binaries with periods longer than about two hours, the loss of angular momentum is 
explained by standard astrophysical theory (interaction between the magnetic field and the stellar 
wind of the secondary). But for short-period binaries, with period less than about two hours, the 
only way to explain the loss of angular momentum is to invoke gravitational radiation. Now, from 
the quadrupole formula, we have 

32G^ M1M2M 
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Inserting this into (§) one can then predict what should be M2 in order that d/a ~ 0, and the result 
is in good agreement with the mass transfer measured from the X-rays observations of cataclysmic 
binaries. 

Another important aspect of the theory of gravitational radiation, with obvious implications in 
astronomy, is the interaction of the gravitational wave field with photons. In this article we present 
the results of a thorough investigation of light propagation in the gravitational wave field generated 
by some isolated system. Our motivation is that electromagnetic waves are still the main carrier of 
astrophysically important information from very remote domains of our universe. Also, the operation 
of interferometric gravitational wave detectors and other techniques used for making experiments in 
gravitational physics (lunar laser ranging, very long baseline interferometry, pulsar timing, Doppler 
tracking, etc.) are fully based on the degree of our understanding of how light propagates in variable, 
time-dependent gravitational fields generated by various celestial bodies. Although quite a lot of 
work has been done on this subject (see, for example, ||2^-|Q) a real progress and much deeper 
insight into the nature of the problem has been achieved only recently 37]-[p9[. The main advantage 



of the integration technique which has been developed for finding the light-ray trajectory perturbed 
by the gravitational field is its account for the important physical property of gravitational radiation. 



namely, its retardation character. Previous authors, apart from Damour and Esposito-Farese |40|, 
accounted for the retardation of the gravitational field only in form of plane gravitational waves. 
Hence, effects produced in the near and induction zones of isolated astronomical sources emitting 
waves could not be treated in full detail. As a particular example of importance of such effects we 
note the problem of detection of gravitational waves created by g-modes of the Sun. The space 
interferometer LISA will be able to detect those waves. However, the problem is that LISA will 
fly in the induction zone of the emission process of these gravitational waves and, hence, a much 
more complete theoretical analysis of the working of the detector is needed. The approximation of 
a plane gravitational wave for the description of the detection procedure is definitely not sufficient. 
The other example could be effects caused by the time-dependent gravitational field of the ensemble 
of binary stars in our galaxy. Timing of high-stable millisecond pulsars might be a tool for the 



detection of stochastic effects produced by that field |41| 



2 Wave generation from isolated sources 
2.1 Einstein field equations 

The gravitational field is described in general relativity solely by the metric tensor g^j^ (and its 
inverse g^'^)- It is generated by the stress-energy tensor of the matter fields T^^ via the second- 



order differential equations 

R^^'^ _ -g^'''R = ^ T^^ , (4) 
2 c 

wliere R^^^ and R = gpaR^" denote, respectively, the Ricci tensor and the Ricci scalar. We assume 
that the matter tensor T^'^ corresponds to an isolated source, i.e. T^'^ has a spatially compact 
support with maximal radius a, and that the internal gravity of the source is weak in the sense that its 
mass M satisfies GM <C ac^. Within these conditions it is appropriate to write the metric g'^'^ in the 
form of a small deformation of the flat metric r]^" = diag(— 1, 1,1,1). We pose h^^ = \/—gg^'^ — r}^'^ 
{g = determinant of g^^y) and assume that each component of h^'^ is numerically small: \h^^^\ <C 1. 
We lower and raise all indices of our metric perturbation h^^'^ with the flat metric; for instance, 
h^v = Vt^pViycrh^'^ ™d h = rjP^hptj. Then the field equations (§) can be re- written in terms of the 
metric perturbation h^'^ by separating out a second-order linear operator acting on h^'^ ^ and the 
remaining part of the equations, which is at least quadratic in h^^ and its first and second derivatives, 
we conventionally set to the right side of the equations together with the matter tensor. This yields 

□ /i^'^ - d^H'' - d'^B^ + ri'^'dpHP = r^'^ , (5) 

where □ = Ojj denotes the flat d'Alembertian operator and where = d^h^'^; on the right side of 
the equation we have put 

r^- = {-g) T^^- + -A^ AM- , (6) 
IdttG 

which represents the total stress-energy distribution of both the matter flelds - flrst term in (P) - and 
the gravitational fleld itself - second term involving the non-linear gravitational source A'^'^ = 0{h?') 
(note that t^^ transforms as a Minkowskian tensor under Lorentz transformations). The divergence 
of the left side of @ is identically zero by virtue of the Bianchi identities, therefore the pseudo-tensor 
T^^ is conserved in the ordinary sense, 

duT^" = , (7) 

which is equivalent to the covariant conservation of the matter tensor, V^T^'^ = 0. A gauge 
transformation h^'^ — > h^^'^+d^^^+d'^£^^—rj^^'^d\^^ does not affect the left side of (|5|), and consequently 
by solving for a vector .^^ the wave equation = —H^^ one can arrange that h^^ satisfies the 
harmonic- gauge condition duh^'^ = 0. In this gauge the field equations (|5|) simply become 

□ h>"' = T^"" . (8) 

We want now to formulate the condition that the source is really isolated, i.e. it does not 
receive any radiation from other sources located far away, at infinity. Recall that we can express 
any homogeneous regular solution of the wave equation n/ihom = at a given field point in terms 
of the values of /ihom at some source points forming a surrounding surface at retarded times. This 
is the Kirchhoff formula (see e.g. [^]), which reads in the case where the surrounding surface is a 
sphere, 

47r 

where p = |x — x'| and t = t' — p/c. To formulate the no-incoming radiation condition we say that 
there should be no such homogeneous regular solutions /ihom (since they correspond physically to 
waves propagating from sources at infinity). Taking the limit r +oo with t + r/c= const, in 
Kirchhoff's formula, we then arrive at the physical conditions that 



/ihom(x',t') = j j 



-Q^KPhhom) + -^[Phhom] 



(x,t), (9) 



lim 



r — '■ + 00 

t+-=const. 



(x,*)=0, (10) 



and that rd\h^^^ should be bounded in this limit. The no-incoming radiation condition (|To|) is thus 
imposed at (Minkowskian) past null infinity J~ in a conformally rescaled space-time diagram. 



2.2 Multipole expansion in linearized gravity 



For the rest of this Section (and also in Section 5) we shah restrict ourselves to the case of linearized 
gravity, defined in particular by the neglect of the non- linear gravitational source term A^*^, for 
which the field equations in harmonic gauge di^h^^'^ = read 

□ h^"' = r^'^ . (11) 

cr 

Within the linearized approximation the matter stress-energy tensor is divergenceless: d^T^^ = 0. 
Therefore the linearized approximation is inconsistent as regards the motion of the matter source, 
which in this approximation stays unaffected by the gravitational field. However this approximation 
is quite adequate for describing the generation of waves by a given source (for instance acted on by 
non-gravitational forces). Prom the no- incoming radiation condition (|lO|), we find that the unique 
solution of (|l^) is the retarded one: 

/.^'^(x,t) = --^y ^— -^r^-(x',t--|x-x'|). (12) 

Since we are being interested in the wave-generation problem, we choose the field point outside 
the source, that is r = |x| > a (with the origin of the spatial coordinates at the center of the ball 
with radius a, so a > |x'|), and we decompose in that region /i^*^ into "multipole moments". The 
straightforward way to do this is to employ the standard Taylor formula for the formal limit x' — > 0, 



1=0 



Ti^',t-r/c) 



(13) 



Notice the short-hand notation L = iii2 ■ ■ ■ ii for a multi- index with / indices, as well as = x'^ = 
^/n^/«2 . . . x'^' , Ol = di-^di^ ' ' ' where di = d/dx^. From this Taylor expansion we immediately 
arrive at the following expression for the multipole decomposition of the metric perturbation. 



„ +00 ,_y 
1=0 



-wr(t--) 
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(14) 



where the "multipole moments" depend on the retarded time u = t — r / c and are given by 

Wf(n) = j (iVx^r^'^(x',u), (15) 



In ( |T^ we employ the notation M. to distinguish the multipole expansion M.{h) from h itself. Of 
course, we have numerically M{h) = h outside the source, however inside the source Ai{h) and h 
will differ from each other; indeed /i is a smooth solution of the equations (Il2|) while A4{h) satisfies 



□A^(/i) = and becomes infinite when r — > [as it is clear from (14)]. In Section 4, dealing with 
the general case of the non-linear theory, we shall prefer to use the multipole expansion in terms 
of symmetric and trace-free (STF) multipole moments. In the present case it is simpler to use the 
non-STF moments Ti'^ . A systematic investigation of the STF multipole expansion in linearized 
gravity can be found in ||43|| . 

Applying the (linearized) conservation law dyT^^ = we easily find certain physical evolution 
equations (and conservation laws) to be satisfied by the multipole moments (^). Making use of 
the Gauss theorem to discard spatial divergences of compact-support terms we successively obtain 
(with (n) referring to n time-derivatives) 



1 



(1) 



-nf = m^;%, (16) 



1 

-^nf = /(/-i)wtY, (1^) 



where the round brackets around spatial indices denote the symmetrization (and where L — 1 = 
ii • • • ii-i] L — 2 = ii - ■ ■ ii-2)- As a consequence of ( [l6| ) we see that the anti-symmetric part of W-^ 
in the indices ij is constant. A more general consequence is 

In the case of the lowest-order (/ = and / = 1) multipole moments the right sides vanish, and 
therefore these equations represent the conservation laws for the corresponding moments. Over 
all we find ten conservation laws, one for the mass-type monopole or total mass M, three for the 
mass-type dipole or center of mass position Xi (times M), three for the time derivative of the mass 
dipole or linear momentum Pj, and three for the current-type dipole or total angular momentum 
Si- Specifically, we define 
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-^/t — 




(19) 


Pi = 
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MXi = 


1 -1/00 _ 

2 * 


r rpOO 
1 d ^ Xi . 
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(22) 



Then, from (|T^)-(|T^, we have 

M = 0, Pi = 0, Si = 0, and Xi = ^ . (23) 



3 The quadrupole moment formalism 
3.1 Multipole expansion in the far region 

We analyze the gravitational field in the far-zone of the source, in which we perform the expansion 
of the multipolar expansion M(h^'^) when r +oo with t — r/c =const (Minkowskian future null 
infinity J^~^). To leading order 1/r the formula (|l^ yields 



^(^n = -^E57i^rN + or^ , (24) 

1=0 ' ^ ^ 



[where (l) represents the Ith time-derivative and where = = n^^n^'^ • • • n*' with n* = x^/r]. 
Because of the powers of 1/c in front of each multipolar piece, it is clear that the far-zone expansion 
of A4 (h) is especially useful when the numerical value of each term of the formula (^) really scales 
with the factor 1/c' in front of it. This will be the case when the typical velocities of the particles 
composing the system are small with respect to the speed of light, z;/c = e<Cl, or, equivalently, 
when the maximal radius a of the system is much smaller than the wavelength A of the emitted 
gravitational radiation (A = cP where P is the typical period of the internal motion). In particular, 
this "slow motion" assumption is always realized in the case of a self-gravitating system with weak 
internal gravity, for which we have 



V 



a GM 



r\j — r\j 



c A V c^a 



<. 1 . (25) 



Thus, for slowly-moving systems we can retain only the first few terms in the multipolar-post- 
Newtonian expansion (p3). Let us restrict ourselves to the terms 
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(26) 

(27) 
(28) 



For easier notation we do not indicate the multipole expansion A4, nor the neglected 0(l/r^) in 
the distance to the source. Using the conservation laws (|T^)-([T7|) we can easily re-express the latter 
expressions in terms of the total mass M, total linear-momentum Pj, and the Newtonian quadrupole 
moment 



Qij 



r2 



00 



2^00 

—^XiXj 



(29) 



Since M and Pi are conserved they do not participate to the radiation field which is therefore 
dominantly quadrupolar. Restoring the neglected post-Newtonian error terms we obtain 



l^'. + lL (n) + 0(e^)}, 



c^r\2 ^ 



(n) + 0(e) 



(30) 
(31) 
(32) 



Note that the contribution of the angular momentum Si appears only at the sub-dominant order 
0(l/r^) [see for instance ( [78[ ) below]. When acting on terms of order 1/r in the distance like in 
(|30|)-(|32|) the derivative di, is proportional to the (Minkowskian) null vector k,^ = (— l,n); namely 
du = —kudo + 0(l/r^). Using this, the harmonic gauge condition d^h^'^ = reads 



(1) 



ku h>"'= O ( ^ 



1 



(33) 



which is checked directly to be satisfied by the expressions (pO|)-(p^). 
3.2 The far-field quadrupole formula 

Using the gauge freedom h'^'^ — > /i^'^ -|- d'^^'^ + d^S^^ — r]^^d\^^, we apply a gauge transformation to 
what is called the transverse-traceless (TT) gauge. Namely we pose 



e = 
e = 
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(2) (2) (2) ■ 

riiab Qab +^1 Qaa -^'^a Q ia 



(34) 
(35) 



The new metric in TT coordinates, say /iJJ (where we lower indices with the flat metric), is straight- 
forwardly seen to involve in its 00 and Oi components only the static contributions of the mass 
monopole and dipoles, namely 



(36) 



In the TT gauge the only radiating components of the field are the spatial ones, ij, and we obtain 



f'lj^ = -^PiMn)l^lb (n) + O + O (^1) . (37) 

The latter equation is known as the far-field quadrupole equation. The TT projection operator 
is defined by Pijab = PiaPjb — ^PijPab, where Pij = 6ij — niUj. This is a projector: namely 
PijkiPklab = Pijabi onto the plane orthogonal to n: thus for instance, riiPijab = 0; and it is trace free: 
P'ijab^ab = 0, SO wc have substituted in ( |37D the trace free part of the quadrupole moment Qab-, i-e. 



hj = Q^j -^Qaa+O (e^) . (38) 

We have added a remainder O(e^) to indicate the post-Newtonian corrections in the source moment 
II computed in Section 4 [note that the remainder in (|37|) is only 0(e)]. 

As we see all the physical properties of the gravitational wave are contained into the TT projec- 
tion (^7|). As a consequence, the effects of the wave on matter fields are transverse: the motion of 
matter induced by the wave takes place only in the plane orthogonal to the propagation of the wave. 
Furthermore, from the trace-less character of the wave we see there can be only two independent 
components or polarization states. We introduce two polarization vectors, p and q, in the plane 
orthogonal to the direction of propagation n, forming an orthonormal right-handed triad. In terms 
of these polarization vectors the projector onto the transverse plane reads Pij = pipj + QiQj- The 
two polarization states (customarily referred to as the "plus" and "cross" polarizations) are defined 
by 



, PiPj ~ QiQj ,TT I PiQj ^ PjQi ,TT /onN 

h+ = — hj , hy, = — ^ hj . (39) 



Until very recently all expectations was that any astrophysical (slowly- 
moving) source would emit gravitational radiation according to (at least dominantly) the quadrupole 



formula (37), involving the mass-type quadrupole moment lij. For instance the waves from the 
binary pulsar obey this formula. However, it has been realized by Andersson and Friedman and 
Morsink [45| that in the case of the secular instability of the r-modes (rotation, or Rossby modes) of 



isolated newly-born neutron stars, the gravitational radiation is dominated by the variation of the 
current-type quadrupole moment Jij. Here we give, without proof, the formula analogous to (^) 
but for the current quadrupole: 

^Jfcurrent = ^ Pijabirv) eacdUc Jbd {u) , (40) 

where the (trace-free) current quadrupole moment is given by 

d^^ Xj)Xa— + O {e'^) . (41) 

It could be, rather ironically, that the first direct detection of gravitational waves would follow the 
formula (Hol) rather than the classic formula (p7|) appearing in all text-books such as P]. 



3.3 Energy balance equation and radiation reaction 

The stress-energy pseudo-tensor of all matter and gravitational fields (in harmonic coordinates) t'^'^ 
is defined by Now, for gravitational waves propagating in vacuum at large distances from their 
sources (in regions where the waves are almost planar), it is appropriate to define the stress-energy 
tensor of the waves as the gravitational source term [involving A^'^(h)] in the definition of t^'^, in 
which hfj^i, is replaced by the far- field metric (|37[). Since the expression of the metric is valid up to 



fractional terms 0(l/r^) in the distance, and since A^'^ is at least quadratic in h, the stress-energy 
tensor of gravitational waves will be valid up to 0(l/r^). Thus, we define, in the far-zone. 

Now to quadratic order A.^^'^ is a complicated sum of terms like hddh + dhdh. But when using 
(33) together with the fact that /c^ = 0, this sum simplifies drastically and we end up with [still 



neglecting 0(l/r^)] 



2 (1) (1) 2 

T^"^ - _2. k^'k'' h'^'^h'^'^ - — A-'^P 



(1) , (1) 



(43) 



The second form is obtained from the inverse of (^9|): hj^ = (piPj — qiqj)h+ + {pigj +Pjqi)hx- The 



expression ( |43| ) takes the classic form a k^k^ of the stress-energy tensor for a swarm of massless 
particles (gravitons) moving with the speed of light. Notice from ( ^3[ ) that the energy density of 
waves is positive definite. In the general case where we do not neglect the terms 0(l/r^) the previous 
expressions of are still valid, but provided that one performs a suitable average over several 
gravitational wavelengths (see ||8|). For quadrupole waves, substituting the quadrupole formula (^), 
we get 



G (3) (3) 

^GW = k'^k'^nM^) I / ki . (44) 



We can integrate the conservation law d^r^" = over the usual three-dimensional space (volume 
element d?x), and use the Gauss theorem to obtain a flux of through a surface at infinity 
(exterior surface element dS'i), so that 



d^x = -c / dSi T^ir ■ (45) 



d 
di 

We consider the n = component of this law, substitute for the expression ( ^^ at the 

quadrupole approximation, perform the angular integration assuming for simplicity a coordinate 
sphere at infinity (i.e. dSi = r'^nidQ), and obtain the famous Einstein (mass- type) quadrupole 
formula 



dE _ G { (3) (3) 



|/^,/., +0(e2)| , (46) 



where E = J d^x t^^ represents the energy (matter + gravitation) of the source, and where we 
re-installed the correct post-Newtonian remainder O(e^). Without proof we give also the formula 
corresponding to the current-type quadrupole moment, 



dE\ 



current 



16G (3) (3) 

'45^ "^^J-^^i ' 



(47) 



where the current quadrupole moment is defined by (41). 

Interestingly, we can treat the decrease of the energy as the result of the back-action of a radiative 
force (cf. [|B]). We operate by parts the time-derivatives in (H^) so as to obtain 



5^5 



G (1) (5) 
/ ij I ij 



5c5 



(48) 



where we put on the left side a term in the form of a total time-derivative, representing a correction 
of order to the energy E, given by 



6E. 



G_ 

5c5 



(3) (2) 
I ij I ij 



(4) (1) ■ 
/ ij I ij 



(49) 



Now, after a time much longer than the characteristic period of the source (for definiteness one 
can consider a quasi-periodic source or perform a suitable average; see e.g. [0), the contribution 
due to the correcting term (|4^) will become negligible as compared to the right hand side of (^). 
Therefore, in the long term, we can ignore this term and finally, the equation ( ^8| ) can be re-written 
equivalently in a form where the energy loss in the source is the result of the work of a radiation 
reaction force Freao namely 



dE 
'dt 



d^X F,eac • V (50) 



where 



F:,,,it, ^) = — px, I (t) , p ^ roo/c^ . (51) 

The equation (^0|)-(^) is called the radiation-reaction quadrupole formula; the specific expression 
(51) of the radiation reaction force is called after Burke and Thorne ^, This force is to 



be interpreted as a small Newtonian-like force superposed to the usual Newtonian force at the 
2.5PN order (or e^). Actually, the Burke-Thorne radiation reaction force is valid only in a special 
gauge. That is, only in a special gauge, differing for instance from the harmonic or ADM gauges, 
does the source equation of motion involve at 2.5PN order the correcting force (pC|)-(|5l]). (See [48| 



for a discussion of various expressions of the radiation-reaction force in different gauges.) Notice 
that the reaction force (|5l| ) contains time derivatives up to the fifth order inclusively. In practice, 
for implementation in numerical codes, high time-derivatives have the tendency of decreasing the 
precision of a numerical computation, and therefore it is advantageous to choose other expressions 



of the reaction force for implementation in numerical codes ||4^, 5C]. On the other hand, the order of 
second, third, and higher time-derivatives can be reduced by making use of the Newtonian equations 
of motion of the matter source. Subsequent implication of such a form of the radiation reaction in 
binary systems leads, for example, to the theoretical prediction of the rate of orbital decay shown 
in (I). 

4 Post-Newtonian gravitational radiation 

4.1 The multipole moments in the post-Newtonian approximation 

In Section 2 we presented the formula for the multipole expansion of the field outside the source in 
linearized gravity. In the present section let us present, without proof, the corresponding formula 
in the full non-linear theory, i.e. when the Einstein field equations (^ are solved taking into 
account the gravitational source term K^^^ . The formula will be valid whenever the post-Newtonian 



expansion is valid, i.e. when (25) holds. Under this assumption the field in the near-zone of a 



slowly-moving source can be expanded in non-analytic (involving logarithms) series of 1/c |1?]. The 
general structure of the expansion is 

F'^(t,x,c) = 5]^/i^,^(f,x), (52) 



9,1 



where hpq are the functional coefficients of the expansion (p, q = integers, including the zero). The 
general multipole expansion of the metric field M{h) is found by requiring that when re-developed 
in the near-zone in the limit of the parameter r/c — > (which is equivalent with the formal re- 
expansion in the limit c oo) it matches with the previous post-Newtonian expansion ( ^2|) in the 
sense of the mathematical techniques of matched asymptotic expansions, i.e. 



M{h)=M{h) 



(53) 



It is worthwhile noting that the equaUty (^) should be true in the sense of formal series, i.e. term 
by term in each coefficient after both sides of the equation are re-arranged with respect to the same 
expansion parameter. 



We find [21, 51] that the multipole expansion generalizing (|lj) to the full theory is composed of 
two terms, 

Mih^^n = finite part a],'[M{An] " — E 7^'"^* " "/"^ ' ^^^^ 

where is the inverse fiat space-time retarded operator. Herein, the first term is a particular 
solution of the Einstein field equations outside the matter compact support, i.e. it satisfies D/ipart — 
A^'^, and the second term is a retarded solution of the source-free (homogeneous) wave equation, i.e. 
'— '^hom ~ "^^^ "multipole moments" parametrizing this homogeneous solution are given explicitly 
by an expression similar to (15), 



H'^^'iu) = finite part J d^x XLTf"'{x,u) , (55) 

but involving in place of the matter stress-energy tensor T^'^ the post- Newtonian expansion t^'^ , in 
the sense of (^2|), of the total (matter+gravitation) pseudo-tensor t^^'^ defined by @. Both terms 
in ( |5^ ) involve an operation of taking the finite part. This finite part can be defined precisely by 
means of an analytic continuation (see for details), but in fact it is basically equivalent to taking 
the finite part of a divergent integral in the sense of Hadamard [^]. Notice in particular that the 
finite part in the expression of the multipole moments ( ^5|) deals with the behaviour of the integral 
at infinity: r ^ oo (without the finite part the integral would be divergent because of the factor 
XL ~ in the integrand and the fact that the pseudo-tensor r^*^ is not of compact support). One 
can show that the multipole expansion (|54[)-(^) is equivalent with a different one proposed recently 



by Will and Wiseman |p2| 



Generally, it is more useful (for applications) to express the multipole expansion not in terms 
of the moments (|55|), but in terms of symmetric trace- free (STF) moments. We denote the STF 
projection with a hat, xl = STF(x^), so that, for instance, Xij = XiXj — |(5jjX^. Then it can be 
shown that the STF multipole expansion equivalent to (^)-(|55|) reads, 



M{hn = finite part □^i[A^(A^'^)] - ^ E r"^^'^^* " "^'""^ ' ^^^^ 



1=0 

where the parametrizing multipole moments are a bit more complicated, 



; j d^x XL j 



T^(u) = finite part / d'^x xl / dz 6i{z)t'"'{x, u + z\x\/c) . (57) 



With respect to the non-tracefree expression (^5]) this involves an extra integration over the variable 
z, with weighting function 



(2^ + l)'!n .2,1 



1 



H^)= \iiif {^-^'y, J_^dz6i{z) = l, hmjiiz) = 6{z) . (58) 



The results (56)-(p8|) permit us to define a very convenient notion of the source multipole mo- 
ments. Quite naturally, these are constructed from the ten components of J-j^'^{u). First of all, we 
reduce the number of independent components to only six by using the four relations given by the 



harmonic gauge condition dyh^^^ = 0. Next we apply standard STF techniques (see 43, 51| for 
details), and, in this way, we are able to define six STF-irreducible multipole moments, denoted 
II, Jl-,Wl,Xl,Yl, Zl, which are given by explicit integrals extending over the post-Newtonian- 
expanded pseudo-tensor t^^ like in (^) . All of the moments Il,Jl,-'' iZl are referred to as the 



moments of the source; however notice that, among them, only the moments (mass- type mo- 
ment) and Jl (current-type) play a physical role at the linearized level. The other four moments 
Wl,Xi, Yi, Zi simply parametrize a linear gauge transformation and can often be omitted from the 
consideration. Only at the order 2.5PN or do they start playing a physical role. The complete 
formulas for the moments Jl are |51] 

Il{u) = finite party d^x j dzlSixiT. - ^^^^ ^ ^^^^^ ^ ^^ Si+iXiLdtT,i 



2(2/ + 1) 

Jl{u 









1 dz EahKii 1 







2/ + 1 

6i+iXL-i>acdt^bc}{'^,u + z\:k.\/c) . (60) 



c2(/ + 2)(2/ + 3) 

In these expressions, <> refers to the STF projection, and we have posed 



7^00 _j_ 7j=0i 



C2 



(where = Si,r^) , J^i = — , E^, = . (61) 



The moments II, Jl given by these formulas are valid formally up to any post-Newtonian order. 
They constitute the generalization in the non-linear theory of the Newtonian moments introduced 



earlier in (|38D and (|4l|). 

In order to apply usefully these moments to a given problem, one must find the explicit expres- 
sions of the moments at a given post-Newtonian order by inserting into them the components of 
the pseudo-tensor T^'^ obtained from an explicit post-Newtonian algorithm. Without entering into 
details, we find for instance that at the IPN order the mass-type source moment II is given (rather 
remarquably) by a simple compact-support formula |17, 21], on which we can, therefore, remove the 
finite part prescription: 

T /"^3 / - , IxP^L ^2 4(2/ + l)xiL „ 1 , / 4^ /«o\ 



We denote the compact-support parts of the source scalar and vector densities in (|61| ) by 

2^00 _j_ rpii q^Oi 

a = 5 , CTj = . (63) 



See Blanchet and Schafer |53[ for application of the formula (|62| ) to the computation of the relativistic 
correction in the Pf, of a binary pulsar [given to lowest order by (|l|)]. On the other hand, Damour, 



Soffel, and Xu ||5J] (extending previous work of Brumberg and Kopeikin |55], |£6[) used the formula 
in their study of the Solar-system dynamics at IPN order. The property of being of compact 
support is a special feature of the IPN mass-moment 1^. To higher-order (2PN and higher) the 



mass moment II is intrinsically of non-compact support (see its expression in |21]); hence the finite 



part prescription in the definition of the moment II plays a crucial role at 2PN. Similarly, starting 
already at IPN order, the current-moment Jl is intrinsically of non-compact support [p!9| , [2l| ] . 

In a linear theory, the source multipole moments coincide evidently with the radiative multipole 
moments, defined as the coefficients of the multipole expansion of the 1/r term in the distance to 
the source at retarded times t — r/c =const. (this is evident from Section 2). However, in a non- 
linear theory like general relativity, the source multipole moments interact with each other in the 
exterior field through the non-linearities. This is clear from the presence of the first term in (^), 
containing the gravitational source A^'^ , and which does contribute to the 1/r part of the metric at 
infinity. Therefore the source multipole moments must be related to the radiative ones, the latter 
constituting in this approach the actual observables of the field at infinity. 



In the TT projection of the metric field one can define two sets of radiative moments (mass- 
type) and Vl (current-type). The definition of these moments is that they parametrize the 1/r-term 
of the ij components of the metric in TT gauge. Thus, extending the formula (|37|), the radiative 
moments are given by the decomposition 



AG 



l>2 



nL-2UabL~ 



21 



b)dL-2 



(64) 



The radiative moments Ul,Vl are related to the l-th time-derivatives of the corresponding source 
moments. Let us give, without proof, the result of the connection of the radiative moments to 
the source moments (p9|)-(|60|) to the order (or 1.5PN) inclusively. To this order some non-linear 
"monopole-radiative ^-pole" interactions appear, which correspond physically to the scattering of the 
l-pole wave on the static curvature induced by the total mass of the source (i.e. the mass monopole 
I) - an effect well known under the name of tail of gravitational waves. We find [pH, 21 1 
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(65) 
(66) 



The same expressions come out from the Will and Wiseman formalism |22| . Here, 6 is a normalization 
constant (essentially irrelevant since it corresponds to a choice of the origin of time in the far zone), 
and K;, TT; are given by 

2^2 + 5/ + 4 ^1 l-l ^1 



l{l + l){l + 2) 



(67) 



k=l " ' k=l 

The first term in ( |65| ) relates essentially to the original quadrupole formula (see (|37|)). From (|65|)-(|6f 
one sees that the first non-linearity in the propagation of the waves is at 1.5PN order with respect 
to the quadrupole formula. Thus, with enough precision, one can replace in ( |6^ ) the mass-type 
moment II by its compact-support expression that we computed in (^) [since the post-Newtonian 
remainder in (|6^) is 0(e^)]. 



4.2 Post-Newtonian radiation reaction 

Emission of gravitational radiation affects the equations of motion of an isolated system dominantly 
at the 2.5PN order beyond the Newtonian acceleration. In a suitable gauge the radiation-reaction 



force density at the 2.5PN order is given by the quadrupole formula (51). In this Section we extend 
this quadrupole formula to include the relativistic corrections up to the relative 1.5PN order, which 
means the absolute 4PN order with respect to the Newtonian force. The method is to compute the 
radiation reaction by means of the matching [in the sense of ([5^)] of the post-Newtonian field to the 
exterior multipolar field. Indeed, recall that the post-Newtonian field is valid only in the near zone, 
and, thus, only via a matching can it incorporate information from the correct boundary condition, 
viz the no- incoming radiation condition imposed at infinity by equation ([lO|), which specifies the 
braking character of gravitational radiation reaction. 

To the relative 1.5PN order, and in a suitable gauge, it can be shown that the reaction force 
derives from some "electromagnetic-like" scalar and vector reaction potentials y^^^^ and y^^^^. 
Explicitly we have |57| 



y-(x,t) = 



(5) , ^ 4GM 
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(68) 



The dominant term in the formula (|6^)-(|69|) is the standard Burke-Thorne reactive scalar potential 
at 2.5PN order [compare ( |68[ ) with ([5l|)]. In this term, consistently with the approximation, one 
must insert the IPN expression of the moment as given by (^). The Burke-Thorne term is of 
"odd" -parity- type as it corresponds to an odd power of 1/c, and thus changes sign upon a time 
reversal (or more precisely when we replace the retarded potentials by advanced ones). Similarly 
is the next term in (|68[)-(|69|) is 3.5PN (i.e. e^), involving both the mass-quadrupole, mass-octupole 
and current-quadrupole moments (the term in the vector potential V^^'^ corresponds really to 

in the equations of motion). However, notice that the next term in the reaction scalar potential 
y^*^*^, at 4PN or order, belongs to the "even" -parity- type. Nevertheless this term is really part of 
the radiation reaction for it is not invariant under a time reversal, as is involves an integration over 
the "past" history of the source, so that when changing the retarded potentials to advanced ones 
the integration range would change to the whole "future"; hence, this term does not stay invariant. 
It represents the contribution of tails in the radiation-reaction force, and is nicely consistent, in the 
sense of energy conservation, with the tails in the far zone [equation (|65|)]. For explicit computations 
of the back-reaction to 3.5PN order in the case of point-mass binary systems see Iyer and Will ||5^, ^ , 
and Jaranowski and Schafer pO| . 

Using the matching (^) one finds that the near-zone post-Newtonian metric (to 1.5PN relative 
order in both the "damping" and "conservative" effects) is parametrized in this gauge by some 
generalized potentials 



reac 



(70) 



The first term represents, to this post-Newtonian order, the conservative part of the metric; it is 
of the normal "even" -parity- type and is given by the usual symmetric integral (half-retarded plus 
half-advanced) of the source densities = (cr, cjj) given by (|6^). The second term denotes 
the radiation-reaction potentials (|6^)-(|6^. By inserting the metric parametrized by (^) into the 
equations of motion of the source (i.e. d^r^'^ = 4^ V^T^'^ = 0), and considering the integral of 
energy, we obtain the balance equation 



dE 

Itt 



(71) 



Here E denotes the energy of the source at the IPN (or even 1.5PN) order. Actually what we obtain 
is not E but some E + SE^/e' + 5Ej/c^ like in ([48|). Arguing as before we neglect these 6E^ and 
SEi. Substituting now the expressions (|68|)-(|69|) for the reactive potentials (and neglecting other 
5E's) we get 
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The right-side is exactly in agreement with the computation of the total flux energy emitted in 
gravitational waves at inflnity, which is computed making use of the stress-energy tensor of grav- 
itational waves (^3|). In particular we recover in the brackets of the first term of (|7^ ) the third 



time-derivative of the radiative moment Uij including its tail contribution. The difference with the 
standard derivation of the flux is that instead of computing a surface integral at inflnity we have 
performed the computation completely within the source, using the local source equations of motion. 



5 Light propagation in gravitational fields of isolated sources 



5.1 General solution of the light propagation equation 



We are going now to calculate in linearized approximation the propagation of a light ray in the 
gravitational field of an isolated source at rest showing a mass-monopole, a spin-dipole and a time- 
dependent mass-quadrupole moment. In linearized approximation, the propagation equation for a 
particle (massless or with mass) with space coordinates reads 

x\t) = -goo,i - goi,t - ■^goo,tx^ - gik,tx^ - {goi,k - gok,i) - (73) 

^^/^^ [ n-i ■ —rij ■ 1 'T^'T^ -\- [ —m ■ , nrxi 1 



goo,kX X - ygikj - -gkj,i jx + y-gkj,t - gokj j x x-'x , 

where the metric coefficients g^u = rj^^i, + f^^, in linear approximation, are related with /i^'^ from 
Sections 2-4 through f^^, = —h^^, + ^rj^^h. The dots denote differentiation with respect to time t 
and c = 1 has been put for simplicity. In the linear approximation scheme, the velocity x* appearing 
on the right-hand-side of ( |7^ ) can be treated as a constant vector. For massless particles (photons) 
it has unit length, i.e. = 1. In the following we use = in the right hand side of ([73|). 
The unperturbed motion of photons reads 

x\t) = xi + k\t-to) , (74) 

where Xq denotes the position of the photon at time of emission to. For solving the light propagation 
equation (73) it is very convenient to introduce the new time parameter r defined hy t = t — t*, 



where t* denotes the time of the closest approach of the photon to the source of the gravitational 
field. Then it holds 

x\T)=e + kW, (75) 

where is the vector pointing from the position of the source to the position of the photon at the 
closest approach. Its length is the impact parameter |.^| = d. Obviously, and A;* are orthogonal to 
each other in the euclidean sense, i.e. ^*A:* = 0. Therefore, the length of x*, r = \x\, takes the simple 
form r = + d'^. Introducing the derivatives di = Pijd/d^^ and dr = d/dr, where Pij = 6ij — kikj 
is the projection operator onto the plane orthogonal to fc*, allows the light-propagation equation to 
be written as 

x\t) = h^k''dj^f,-dr(^ky,^ + h^foo-h^k^kPf,p^ , (76) 

where the four-dimensional vector k" reads k" = (1,A;*). 

To get a complete overview of the influence of a gravitational wave, emitted from an isolated 
source, on the propagation of light rays, we use the representation of the metric coefficients (p!4|) 
which is valid all-over in the space outside a domain which includes the matter source. Splitting 
ffj,^, into a canonical part /^^"^ which contains trace-free tensors only, and a gauge part, i.e. ffj,^, = 
+ dfj_Wy + duW^, we obtain in case of mass-monopole, spin-dipole, and mass-quadrupole source 
moments (remember the source being at rest) 
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, (78) 



/r = 5^,/oT + ?4(i-0- (79) 
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Herein, for simplicity, we have put G = 1, and di = d/dx^ . The mass M, spin 5*, and the quadrupole 
moment lij of the source of gravitational waves are given (in the Newtonian approximation) by the 



expressions (|19D, (pi]), and (pq). The explicit expressions for the gauge functions relating f^^^^ 
with ffj^iy are important for general discussion of light propagation in the field of gravitational waves 
emitted by the isolated source. However, for the sake of simplicity they will be omitted. Their 
precise form can be found in [p^ j. 

The insertion of these expressions ([f7|)-([7^) into the equation ( |76| ) results in the equation 

x'{t) = [2M (Bi - kA) - 2SP (^e,p,d,r - A^.^ip,^)] |^| (80) 



ipq ni^pqr ^ ^r^p^iqr j ^ ^ \ 1-3^^ 

where the vector (p^^ denotes terms which are of gauge type. The precise form of ip^ is not important 
here and can be found in [ p7| . 

The solution of equation (|80|) , using the boundary conditions oo) = fe* and x*(to) = Xq - 
emission point in space of the light ray at time tq, reads 



x\t) = k' + E\T), (81) 
x\t) = x%{r) + EHT)-E\ro), (82) 



where x^(t) denotes the unperturbed trajectory ([75D. The relativistic perturbation of the trajectory 
is given by 

H^r) = (2M^i + 2SPkJejp,^ig^A{T,0 + ^^pqBpg{T,0- (83) 
(iMki + 2SPei^ada] 1-1 - f kAa - 2k A ^ ^ ^""^^ ~ 
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(84) 
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whereby the scalar functions A and B and derivatives of the tensors Bij , Cij , Dij are known fully 
explicitly. They are given by 
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P'' + ^) dkBpqir, + P'^diBpqir, + ^dikBpqir, 



(89) 



where the variable y = t — \/t'^ + d'^ is the retarded time argument for the photon which passes 
through the point of closest approach to the source of the gravitational radiation at time t* = 0. More 
details concerning the method of calculation of light ray trajectory in time dependent gravitational 
fields can be found in [37, 



5.2 Time delay and bending of light 

The time delay results in the form 

t-to = |x-xo| 



k • H(r) + k • H(to) 



(90) 



or 



to 



\x - xo\ + AM{t,to) + As{t,to) + AQ{t,to) 



(91) 



where |x — xo| is the usual Euclidean distance between the points of emission, xq, and reception, x, 
of the photon. Am is the classical Shapiro delay produced by the (constant) spherically symmetric 
part of the gravitational field of the deflector (see, e.g. H), As is the Lense-Thirring or Kerr delay 



due to the (constant) spin of the localized source of gravitational waves |39|, and Aq describes an 
additional delay caused by the time dependent quadrupole moment of the source |^. Specifically 
we obtain 
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Let us now denote by a* the dimensionless vector describing the total angle of deflection of the 
light ray measured at the point of observation and calculated with respect to vector given at past 
null infinity. It is defined according to the relationship 



a\T,0 = F[k.H(T,0] 



ir,0 



(96) 



For observers being far away from the source of the gravitational wave the projection of the 
mass-quadrupole tensor of the source of gravitational radiation onto the plane orthogonal to the 
propagation direction of the gravitational wave is the crucial object which enters into the observable 
effects. It reads 
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(97) 



where again we denote ni = x*/r. 

In the case of small impact parameter d {d/rQ <^ l,d/r <C 1) we respectively obtain for the time 
delay and the angle of deflection 



t - to - |x - xol = -Alp + 2M ln(4rro) 



(98) 



where ip is the gravitational lens potential having the form 



^P 
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(Notice that in this gravitational lens approximation n* = fc* holds.) Remarquably, the gravitational 
lens potential does depend on the gravitational source mass-quadrupole tensor only through its value 



I.e. 



at the time of closest approach. Furthermore, the gravitational lens potential decays like 
it is not being influenced by the wave part of the gravitational field |4^, 37|. 

A direct consequence of the time delay formula is the frequency shift formula for a moving 
gravitational source 



6iy ^ dip j 
V at* 
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where is the velocity of the observer. It is worthwhile noting that the expression ( |101D holds 
for the source of electromagnetic waves being at past null infinity. An exhaustive treatment of the 
gravitational frequency shift for arbitrary locations of observer, source of light, and the source of 



gravitational waves is rather complicated and has been done only recently in [38|. 



6 Detection of gravitational waves 

In the asymptotic regime of a gravitational wave field, the time delay reads 



Aq k;t,to = ■ 

1 — cos ( 
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(102) 



where 6 is the angle between the receiver - (light) emitter and receiver - (gravitational wave) source 
direction (cos0 = — iV*A;*) and where the assumption |x — xo| <^ r has been made. Let us now 
apply the above formula to the time delay in a Michelson interferometer. Obviously, in this case, 
r = ro holds. For simplicity we assume that the interferometer device is oriented orthogonal to 
the propagation direction of the gravitational wave, i.e. N^k\ = N^kJ^, where k\ and denote the 
directions of the two interferometer arms which are taken to be orthogonal {klk^ = 0) and of equal 
length L. We also assume that the light, emitted from the beam-splitter, is reflected once at the end 
mirrors. Furthermore, the interferometer arms are to be oriented such that they coincide with the 
main axes of the plus-polarization. Then the relative time delay of the reflected light beams reads 

AQ(ki; t, t-2L) = ^[Ii;^{t - ro) - ij^^it - ro - 2L)] . (103) 

The multiplication of the relative time delay by the angular frequency of the laser light, uj, which is 
treated as constant in the approximation under consideration, results in the measurable phase shift 
at time t of A$(t) = 2u;AQ(ki; t,t — 2L). If we assume for the plus-component of the gravitational 
wave the expression h+{t — ro) = cos{uigt), where cOg is the constant frequency of the wave and 
its constant amplitude, we obtain for the phase shift (cf. ||6T[| ) 

A<^{t) = 2A+— sm{ujgL) cos[ujg{t - L)] . (104) 

LOg 

The maximal amplitude is achieved if the condition ujgL = 7r/2 holds. This yields 

Aina.^>(i) = 2A+— sm{LOgt) . (105) 

UJ„ 



At the photo-diode the following photo-current results |62| 



/ph(t) = /min + ^"'"^ ^ ^'""^ [1 - COS m] , (106) 

where the phase (j){t) is composed out of the signal plus a modulation term from Pockels cells, 
(pmsin{u;mt), i.e. 



(j){t) = A$(t) + (pmSmiuJmt) 



(107) 



The approximate decomposition of the expression (\10(\ ) into a dc ( "direct current" or non-alternate) 
part and a signal part reads /ph(0 = -^dc + lujm^ where 

Idc = /min + ^ [1 - M'Pm) COS A$(t)] , (108) 

Ico^ = IcsJi{(t>m) sin A^>(t) sin(a;„t) , (109) 

with Bessel functions Jq and Ji; Jeff = -^max — -^min- Because of the smallness of the gravitational 
phase shift we get 

Idc = Imin + ^ [1 - MM] , (110) 

= Ic{iJl{(pm)A^it) Sm{uJrnt) . (Ill) 

Using the equation (|105|), the latter equation can be written 



UJg 



lum = IcsJl{(t>m)A+ — { COS [{uJrn " ^g)t\ - COS {{uJm + ^g)t\ \ ■ (112) 



In this side-band form, the signal from the gravitational wave is being detected. For more details 



we refer to [61]. 
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